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Abstract 

We consider the one-dimensional Cattaneo equation for transport of 
scalar fields such as solute concentration and temperature in mass and 
heat transport problems, respectively. Although the Cattaneo equation 
admits a stochastic interpretation - at least in the one-dimensional case 
- negative concentration values can occur in boundary-value problems 
on a finite interval. This phenomenon stems from the probabilistic na¬ 
ture of this model: the stochastic interpretation provides constraints 
on the admissible boundary conditions, as can be deduced from the 
wave formulation here presented. Moreover, as here shown, energetic 
inequalities and the dissipative nature of the equation provide an alter¬ 
native way to derive the same constraints on the boundary conditions 
derived by enforcing positivity. The analysis reported is also extended 
to transport problems in the presence of a biasing velocity field. Sev¬ 
eral general conclusions are drawn from this analysis that could be 
extended to the higher-dimensional case. 

Keywords :Transport phenomena, positivity of operators, hyperbolic equa¬ 
tions, Cattaneo equation, microscopic stochastic models. 

1 Introduction 

Cattaneo transport equation for scalar fields (matter, heat) stands as a 
valuable attempt to overcome the intrinsic problems of Fickian transport, 
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namely its infinite velocity of propagation [T] (throughout this article we 
refer as “Fickian” to the wealth of transport problems where the flux of the 
transported entity is proportional to the concentration gradient, such as in 
the Fick’s equation for mass transport, or to temperature gradient, such as 
in the Fourier equation for conductive heat transfer). 

This is achieved by introducing an exponentially decaying memory con¬ 
tribution in the flux-concentration constitutive equation. The Cattaneo 
model finds important applications in some specific transport problems (such 
as the phenomenon of the second sound in superfluid helium) [2] (see also 
the review by Joseph and Preziosi mm, and in a variety of other, more 
engineering-oriented, studies EEUZj. Some authors have developed an ap¬ 
proximate mesoscopic derivation of Cattaneo-like transport equations in the 
form of a generalized hyperbolic Fokker-Planck equation M, and Cattaneo 
model rests as a fundamental example of generalized constitutive equation 
in the theory of extended irreversible thermodynamics mmm- 

Its implications in theoretical physics should not be underestimated, as 
it represents a valuable contribution for extending transport models in a 
relativistic (Lorentz-invariant) way [Q], and it constitutes the starting point 
for attempting a stochastic interpretation of the Dirac equation mag. 

However, its validity has been deeply questioned and criticized for space 
dimension greater than one, as the associated Green function attains neg¬ 
ative values m, so that the propagation of generic non-negative initial 
conditions can return negative concentration values. 

In point of fact, even this observation can be pushed forward, since there 
exist one-dimensional problems on the finite interval for which the solution of 
the Cattaneo equation in the presence of non-negative initial and boundary 
conditions can become negative. This is shown in Section [2] by means of 
a simple, closed-form example. The Cattaneo equation on the interval has 
been studied also in HUES], for boundary conditions that do not present 
problems as it regards positivity. 

This result may seem in contradiction with the stochastic derivation of 
the Cattaneo model starting from a simple stochastic equation driven by 
dichotomous Poisson noise, developed earlier by S. Goldstein [19] and sub¬ 
sequently re-elaborated by M. Kac in a simple but seminal contribution [20]. 
This contradiction is however purely apparent and can be resolved by for¬ 
mulating the transport problem in the wave-formalism envisaged by Kac 
[20] and subsequently elaborated by several others authors H8l EH [22] |23], 
just to cite some relevant contributions in the field. The wave analysis in¬ 
dicates unambiguously that not all of the boundary conditions, that are 
commonly applied in a Fickian framework, are physically admissible in hy¬ 
perbolic transport schemes. More precisely, the wave-like nature of the 
transport models with memory dictates bounds and constraints on the ad¬ 
missible boundary conditions preserving positivity. This is shown in Sections 
[3] and [4] dealing with the desorption experiment addressed in Section [2] The 
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same constraints derived from positivity emerge from the analysis of the 
L 2 -nornrs, i.e., from energetic conditions, once dissipation is enforced (see 
Section 0. 

In this paper the analysis is also extended to one-dimensional problems 
in the presence of a deterministic biasing velocity held, deriving constraints 
on boundary conditions in the classical problem of boundary-layer polariza¬ 
tion characterizing transport across permeable or perm-selective membranes 
(see Section 0. Some general conclusions, oriented towards a rigorous for¬ 
mulation of hyperbolic transport models, are addressed in the concluding 
Section [71 


2 The Cattaneo equation on the unit interval 


Consider a transport problem for a conserved scalar held u(x, t ) on the 
spatial interval x € [0, L\ and time t, in the absence of source terms. In this 
case, the balance equation is obviously given by 

du(x,t ) dJ(x,t ) 

dt = dx ’ ( ’ 

and assume that the flux J(x, t) is related to u(x, t ) via a constitutive equa¬ 
tion of Cattaneo type 


J(x, t ) = —k 


du(x,t ) dJ(x,t ) 


- Tr 


( 2 ) 


dx c dt ’ 
where k is the “Fickian” diffusivity (conductivity) and r c the characteristic 
memory time. By substituting eq. © into eq. (0 we obtain 


du(x,t ) d 2 u(x,t ) d 2 u(x,t) 

+ = (3) 

namely the Cattaneo equation. 

As t c —> 0, the Cattaneo equation degenerates into a diffusion equation, 
while as t c ^ oo, keeping k/r c = constant, it reduces to a one-dimensional 
wave equation. Both these two limit conditions, i.e., the parabolic limit and 
the pure dissipation-free wave equation, bring physical contradiction inside, 
namely the infinite velocity of propagation characterizing Fickian transport 
on one hand, and the occurrence of pure wave-like propagation without 
dissipation on the other hand, as mentioned by Joseph and Preziosi 121 IH- 
As initial condition at t = 0 let 


u(x, 0) = Ui n = constant, 


du(x, t ) 
dt 


= 0 . 
t =o 


(4) 


Introducing the nondimensional variables y = x/L , 6 = tk/L 2 , a = r c k/L 2 , 
c l i = u/u-ya, the balance equation for cj>(y, 0) reads 


dcj) d 2 4> d 2 4> 
d6 a d0 2 dy 2 ’ 
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and cj)\ e = 1, dc/)/dd\e=o = 0. 

Let us consider a desorption experiment, meaning that u(x,t ) can be 
viewed, for instance, as the concentration of a solute inside some solid matrix 
where it diffuses. Suppose that the boundary at x = 0, (y = 0), is perfectly 
impermeable to transport so that J\ x= o = 0 for all t > 0. Conversely, at 
x = L, (y = 1), solute diffuses out of the solid matrix onto an external 
reservoir that is much larger than the sample matrix and perfectly mixed so 
that the concentration at y = 1 can be assumed vanishingly small, i.e., 

<KM) = 0. (6) 


It follows from the constitutive equation, that the zero-flux boundary con¬ 
dition at y = 0 reduces to 


d<t>(y, 6) 


dy 


= 0 


(7) 


y =o 


From eqs. ©-([71) the concentration field (j>(y,0) can be expanded in Fourier 
series 


e ) = ^ n W cos K.y) > u n = — 9 ^ • 


n= 1 


The Fourier coefficients 4> n (6) satisfy the equation 

Oi <t>n(9) + 4>n(9) + l/ n4 > n{9) = 0 . 


( 8 ) 


(9) 


where <f> n indicates the derivative with respect to 9 and cj) n the second 
derivative, equipped with the initial condition ^> n (0) = 2(— /v n = cj> n fl, 
n = 1,2,... . The characteristic equation associated with eq. © is readily 
A 2 + A /a + = 0, and it admits the solutions 


X = ~'t ± ^ i A n = l-4 v 2 n a. (10) 

Depending on the sign of the discriminants A n two cases occur: (i) for 
n < n*, where A n * > 0 and A n » + i < 0, the roots of the characteristic 
equation are real: Ap„ = —1/2 a + y/A n /2a, A 2 , n = —1/2 a — \J A n /2a. 
Therefore, for n < n*, 

M6) = A n e x ^ e + B n e x ^ d . ( 11 ) 


Enforcing the initial conditions one gets 

A _ ^2,n 4>n ,0 D _ Ai )71 a 

-rt-n — > > } 7J) n — J-\. n . 

^2,n ^1 ,n A2,n 


( 12 ) 


For n > n *, the roots of the characteristic equation are complex conjugate: 
A 12 = —l/2a ± iu n , i = \/—1, u n = y/—A n /2a, and 


<j>n(0) = e e/2a [C n cos (u n 9) + D n sin(w n 6>)] , (13) 
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where: 



( 14 ) 


To sum up, the solution of the above transport problem involving the Cat- 
taneo equation can be expressed as 


n* 



(15) 


n= 1 


oo 


^ [C n cos(uj n 6) + D n sm(u n 9)\ e 9/2a cos(z/ n y) . 


n=n *+1 


Figure Q] depicts the concentration prohles for different values of 9 along 
the nondimensional spatial coordinate y, for a = 0.2 (panel a) and for a = 
0.5 (panel b). In both cases, it can be neatly observed that the concentration 
attains negative values. 

The unphysical inconsistency of the transport model is even more evident 
by considering the relative released fraction 



(16) 


In the case under examination, since the time derivative of the flux at y = 1 
is identically vanishing for 9 > 0, eq. (1161) can be written in the equivalent 
form 





(17) 


similarly to the pure Fickian case. The quantity M d (9)/M 00 represents 
the ratio of solute released up to time 8 to the overall solute quantity 
initially present within the matrix. Since there is no source, physics dic¬ 
tates that this quantity should be bounded by 1, i.e., M d (8)/M 00 < 1, 
linig —>•00 M d (9)/M 00 = 1. 

Figure [2] depicts the time evolution of M d (9)/M Q0 for several values of 
the dimensionless number a, showing the occurrence of an overshot of this 
quantity reaching unphysical values greater that 1. This phenomenon is 
more pronounced as a increases i.e., as the influence of the memory term in 
the constitutive equation becomes more significant. 

It should be observed that overshot dynamics in sorption experiments 
have been observed in transport across polymeric matrices where the vis¬ 
coelastic response of the material is responsible for the occurrence of memory 
effect in mass transport (see [23J and reference therein, and particularly fig¬ 
ure 1 in this article based on closed-form solution of the transmission-line 
equation [25j). However, while in sorption experiments an overshot in the 
overall sorption curve does not violate in principle any fundamental law 
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Figure 1: Concentration profiles cf>(y, 6) vs y for increasing values of 6 = nAO, 
n = 1, 2,..., Ad = 0.1 (in the direction of the arrows). Data are obtained 
from eq. (USD truncating the series up to IV = 40000 modes. Panel (a) refers 
to a = 0.2, panel (b) to a = 0.5. 


of nature, the occurrence of A l,i{9)/M 00 greater than 1 in the desorption 
experiment discussed in this Section is physically inconsistent, since it is 
the macroscopic manifestation of the occurrence of negative concentration 
values. 


6 














































Figure 2: Md(0)/AI oo vs 9 for different values of a. Line (a): a = 1, line 
(b): a = 0.5, line (c): a = 0.2. 

3 Stochastic dynamics, wave formulation and boundary- 
condition constraints 

The paradoxes of negative concentration values and of a released fraction 
greater than one at intermediate timescales provide an interesting example of 
the role played by the boundary conditions, in that the choice of admissible 
boundary conditions is entirely entailed in the field equations, or better to 
say in the present case, in its underlying microscopic stochastic formulation. 
Different classes of boundary conditions for the Cattaneo equation defined 
on an interval have been addressed in mmm- The focus of the present 
analysis is however slightly different with respect to these papers, as it is 
strictly oriented on the consistency of the boundary conditions in order to 
preserve the positivity of the solutions. 

The result obtained in the previous Section could sound as a paradox, 
by considering thta the Cattaneo equation on the one-dimensional line ad¬ 
mits a probabilistic interpretation in term of a stochastic process driven by 
dichotomous Poisson noise [201 HU 1221 [23]. 

Consider the stochastic equation 

dx(t) = c(— l) x ^ dt , (18) 

where c > 0 is a reference velocity, xif) an ordinary Poisson process charac¬ 
terized by the transition rate constant a > 0 [2D]. Let X (t) be the stochastic 
process associated with eq. (fTSl) and let 

p + (x,t)dx = Prob{X(f) € (x, x + dx), (—l) x ^ = 1} 
p~(x,t)dx = Prob{X(f) G (x, x + dx), (— = — 1} , (19) 
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the probabilities of the occurence of a value of X(t) in the interval (x, x+dx), 
and of the stochastic dichotomic perturbation (—1) X W = ±1, respectively, 
at time t. The two partial probability densities p ± (x,t) specify completely 
the statistical characterization of the process (fl8l) [20] . From eqs. (fl8l) - (fT9l) . 
it follows that the evolution equations for the two quantities p±(x,t) read 


dp + (x,t) 

dt 

dp~(x,t ) 
dt 


dp + (x,t ) , . . . 

-c--- ap + (x,t)+ap (x,t) 

ox 

dp~(x, t ) 


dx 


+ ap + (x,t) — ap (x,t) 


( 20 ) 


where a is the rate constant of the Poisson process %(f). By defining the 
two quantities, 


P(x,t) =p + {x,t) +p (x,t), J(x,t) = c\p + {x,t) -p (x,t)] (21) 


it follows that p(x, t ) satisfies the conservation equation 


dp(x,t) dJ(x,t) 

dt dx ’ 

while J(x,t) fulfills the relation 

1 dJ(x,t) T , . c 2 dp(x,t ) 


( 22 ) 


(23) 


The two scalar fields p(x,t ) and J(x,t) correspond respectively to the over¬ 
all probability density function and its flux associated with the stochastic 
dynamics (Hbd, and eqs. (I22l) - ll23l) return the Cattaneo transport equation 
by identifying r c = l/2a and k = c 2 /2a. It follows from the definition, 
that the quantity p(x,t) should be strictly positive due to its probabilistic 
interpretation. Let us refer to p + (x,t ) and p~(x,t) as the partial probabil¬ 
ity densities (or waves), and to the equations (12U1) as the partial probability 
model (acronym PPM). 

Starting from the decomposition into partial probabilities, it is possible 
to develop a theory for the boundary conditions. The partial probabilities 
are the fundamental physical quantities in a transport process on the line 
driven by dichotomous noise, while p(x, t ) and J(x, t ) should be viewed as de¬ 
rived quantities. Consequently, the boundary conditions should be properly 
expressed in terms of p + (x , t) and p~{x , t) to ensure positivity. This concept 
can be expressed in an alternative way: physically admissible boundary con¬ 
ditions for the hyperbolic transport problem on the interval [0, L] are those 
that preserve the positivity of the partial waves p + (x,t), p~(x,t), start¬ 
ing from non-negative (initial and boundary) values of these two quantities. 
The application of this approach to the setting of the boundary conditions 
explains the apparent paradox addressed in Section [2j 










Consider the boundary conditions for the problem addressed in Section 
El namely J\ x= o = 0 and eq. ([6]), i.e., in the present notation p(l,t) = 0. 
The impermeability condition at x = 0 implies for the partial probabilities 

p + {x,t)\ x=Q = p-{x,t)\ x=Q , (24) 

which is a fully admissible boundary condition. Eq. (12411 implies that the 
forward propagating wave (p + ) equals the incoming regressive wave p~ at 
x = 0. Next, consider the other boundary condition eq. JU]). It means that 
[P + 0M) +P~(x,t)] \ X=L = 0, i.e., 

p-(x,t)\ x=L = -p + {x,t)\ x=L , (25) 

implying that, if the forward probability is positive at x = L, the “reflected” 
backward wave should attain negative values. Clearly, it is not an admissible 
boundary condition, the consequence of which is the occurrence of negative 
values for p + (x,t ) +p~(x, t). 

In a similar way, it is possible to address other kinds of boundary con¬ 
ditions. A paradigmatic example is given by radiative boundary conditions 
arising e.g. in interphase transport across a boundary layer. Below, we ana¬ 
lyze this case. The prototype of homogeneous radiative boundary conditions 
is given by 

J{x,t)\ x=L = h p(x,t)\ x=L , (26) 

where h (possessing the physical dimension of a velocity m/s ) is the mass- 
transfer coefficient across a boundary layer at x = L. Substituting the 
expressions for p and J in terms of partial probabilities, eq. (1261) becomes 

P~(x,t)\ x=L = ^ P + (x,t)\ x=L . (27) 

The latter expression indicates that positivity is preserved provided that the 
condition 

c-h> 0 (28) 

is fulfilled. This condition can be expressed in terms of the physical param¬ 
eters k and r c resulting 

De h = ^ < 1 . (29) 

where the dimensionless group Deh can be referred to as the radiative Deb¬ 
orah number. The radiative Deborah number is the ratio of the character¬ 
istic “recombination time” r c (the intrinsic memory timescale) to the dif¬ 
fusion (conduction) time associated with the transfer across the boundary 
layer tu = b 2 /k, where 6 is the boundary-layer width. Since the boundary- 
layer width 5 is implicitly defined by the relation h = k/5 , it follows that 
tu = k/h 2 . 

Eq. (1291) indicates that for the transport problem under consideration, 
radiative boundary conditions are admissible (i.e. they preserve positivity) 
provided that the Deborah number Deh is less than or at most equal to 1. 
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4 Revisiting the boundary conditions 

In the previous Section we have derived the more general form of linear ad¬ 
missible boundary condition in the framework of the one-dimensional PPM: 

p~(y> e )\y=i = 7 p + (y^)\ y= i > ( 3 °) 

where 7 € [0,1] is a dimensionless parameter depending solely on the Deb¬ 
orah number Deh- Enforcing eq. (13011 at y = 1, and the zero-flux boundary 
conditions at y = 0 , it is possible to obtain physically consistent results 
regarding desorption kinetics in hyperbolic transport. 

Figure [3] depicts the overall concentration profile 4>(y, 9) = p(y, 0) = 
p + {y,9) + p~(y,9), solution of the dimensionless PPM, equipped with the 
boundary condition (1301) , 7 = 0 at y = 1 for two different values of the dimen¬ 
sionless parameter a. The initial conditions are p + (y, 0) = p~(y, 0) = 1/2 
uniformly in y € (0,1). As expected from the analysis of admissible bound¬ 
ary conditions developed in Section [3j the occurrence of negative concentra¬ 
tion values is “cured” by eq. {HD- 

The overall desorption kinetics, expressed by the relative release fraction 
M ( j{9) /M^, is depicted in figure [4) Panel (a) refers to 7 = 0 for different 
values of a, while panel (b) depicts the influence of the boundary parame¬ 
ter 7 . As expected, the release kinetics becomes progressively slower as 7 
approaches 1 , since 7=1 corresponds to the impermeability (zero-flux) con¬ 
dition. The data depicted in figure 0] indicate that the unphysical overshot 
of Md{9) /Mqo up to values greater than 1 cannot occur if proper bound¬ 
ary conditions are enforced, accounting for the wave-like propagation of the 
stochastic perturbation. 

With reference to the desorption experiment analyzed throughout this 
work, let 

p(y,9) = l + u(y,9) . (31) 

Since 0 < p(y,9) < 1, it follows that —1 < u(y,9) < 0. In terms of partial 
waves this implies 

P + (y,0) = ^ + u + (y,6) , p~(y,9) = \ + u~(y,9) . (32) 

The balance equations for the partial waves u ± (y,0) are identical to eqs. 
(1201) replacing p± with u , t with 6 and x with y. The initial conditions 
are u + (y, 0 ) = u~(y, 0 ) = 0 , and impermeability at y = 0 is just u + (0,9) = 
u _ ( 0 , 9). 

More interesting is the boundary condition at y = 1. From eqs. m 
and (l32l) one obtains 

u ~(yi 9 )\ y =i =7 u + {y,9)\ y=1 - . (33) 
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Figure 3: Concentration profiles cf>(y , 6) vs y solutions of the PPM with the 
boundary condition (15U1) at y = 1. 7 = 0 for increasing values of 9 = nA9, 
n = 1,2,..., A 0 = 0.1 (in the direction of the arrows). Panel (a) refers to 
a = 0.2, panel (b) to a = 0.5. 


This indicates that the desorption kinetics can be solved by considering the 
complementary partial waves ^(y^O) equipped at the transfer boundary 
(y = 1) with the condition (1331) . This is a non-homogeneous expression 
containing a strictly negative source term —(1 — y )/2 < 0 , since 7 £ ( 0 , 1 ). 

This elementary observation suggests a way to handle boundary condi¬ 
tions in a sorption experiment, which is the complementary problem to the 
desorption kinetics analyzed so far. In a sorption experiment, solute diffuses 
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Figure 4: M^{0)/M oo vs 6 obtained from the PPM. Panel (a): 7 = 0, for 
different values of a. Line (a): a = 1, line (b): a = 0.5, line (c): a = 0.2. 
Panel (b) a = 0.5, for different values of 7 . The arrow indicates increasing 
values of 7 = 0, 0.2, 0.5, 0.8. 


within the solid matrix from y = 1 , where it is in contact with a perfect 
reservoir at constant unit dimensionless concentration. Assume a vanish¬ 
ing initial condition in y € (0,1). The boundary condition at the transfer 
interface y = 1, that is complementary to the desorption condition (1301) is 
simply 

u ~(y^)\y=l =7 u + (y,0)\ y=1 + ■ (34) 
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It corresponds to a change of sign in the source term entering the desorption 
conditions (1331) . Initial conditions are u + | 0 =o = u~\q = q = 0, and at y = 0 is 
enforced. 

In a sorption experiment, the overall uptake kinetics is expressed by the 
sorption curve 



(35) 


Figure [5] (panel a) depicts the concentration profiles for r c = 0.5, 7 = 0.6 
at different time instants sampled at A 8 = 0.1. Panel (b) of the same 
figure compares the resulting sorption curves (solid lines) obtained using 
eq. (1351) and the corresponding desorption curve associated with eq. (T30T) . 
As expected, sorption and desorption curves at the same values of r c and 
7 coincide, confirming that the nonhomogeneous boundary conditions (1341) 
represents in sorption experiments the analogue of the desorption boundary 
conditions m- 

It is important to stress once again that this complementary derives from 
the intrinsic decomposition of the concentration field into partial waves. This 
decomposition, with all of its implications, is the essence of hyperbolic trans¬ 
port. Any attempt in modeling boundary conditions in sorption dynamics 
starting directly from the overall concentration field u(y, 8), thus neglecting 
the underlying decomposition in partial waves 8), leads unavoidably to 

physical inconsistencies. Let us clarify this statement with a simple exam¬ 
ple. In an ideal sorption experiment (starting from u\g = o = 0), the simplest 
boundary condition corresponding to the contact with an ideal reservoir at 
unit concentration would be 


u(y,0) \ y=1 = 1 


(36) 


in the absence of transfer resistances. Theuse of this equation, or of the 
equivalent expression, 



(37) 


ends up with an overshot dynamics for the sorption curve M s (8)/M oa , cor¬ 
responding to local concentration values greater than 1 , i.e. greater than 
the feeding concentration at the boundary y = 1. This is shown in figure O 
depicting the concentration profiles (panel a), and the overall sorption curve 
(panel b) at r c = 0.5. 

5 Norm dynamics and dissipation 

The constraints on the structure of the boundary conditions can be derived 
from the evolution of the integral norms, enforcing dissipation. Consider eq. 
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Figure 5: Panel (a): Concentration profiles in a sorption experiment at 
t c = 0.5, imposing the boundary conditions (PTT1) with 7 = 0.6, sampled at 
Ad = 0.1 (time increases in the direction of the arrow). Panel (b): Sorption 
curves M s (d)/M oa from eq. (|55D (solid lines) for different values of r c and 
7 . Symbols are the corresponding desorption curves Md(d)/M 00 . Line (a) 
and (o): r c = 0.5, 7 = 0.6, line (b) and (•): r c = 0.5, 7 = 0.8, line (c) and 
(□): r c = 0.1 and 7 = 0.6, line (d) and (■): r c = 0.1 and 7 = 0.95. 


( 1201 ) defined on an interval T>. So far no assumption is made for T>, be it 
(— 00 , 00 ) or [0, L\. Multiplying the first equation (|20l) by p + , the second by 
p ~, summing the resulting equations, and integrating over V provides the 
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Figure 6: Panel (a): Concentration profiles u(y,0) in a sorption experiment 
imposing the boundary condition (15S1) . r c = 0.5. The curves are sampled at 
Ad = 0.1 (time increases in the direction of the arrow). Panel (b) Sorption 
curve M s {0)/M oo associated with the data depicted in panel (a). 


expression 


Sv^^ 2 ” (r?] ^ 

a / ( P + — P~) 2 dx . 

Jv 


( 38 ) 
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Let us express the quantities entering eq. (13811 in terms of concentration and 
flux: 

(; p + ) 2 + (p~ f = \{p + +p~f + \{p + -p~) 2 = \p 2 + ^2 j2 
0 p + ) 2 -{p~) 2 = {p + +p~)(p + -p~) = \p j ■ ( 39 ) 

Two cases can be considered: either V = (— 00 , 00 ) or V = [0, L\. 

On the real line, i.e. T> = (— 00 , 00 ), regularity conditions at infinity 
apply, and consequently the first integral at the r.h.s. of eq. ([551) vanishes 
identically. 

Letting ||/||| 2 (i) = fp f 2 ( x > t) dx be the square L 2 -norm of a square- 
summable real-valued function f(x,t) in V , eq. (15H1) can be rewritten as: 

j t (iwii. w + ^ Piii.w) = -\ mis ). (40) 

which expresses the dissipation properties of hyperbolic “diffusion” i.e. of 
the stochastic motion driven by pure dichotomous Poisson noise. In the limit 
c —>• 00 , J —>• —kdp/dx, eq. (HU1) reduces to the classical norm condition, 
d\\p\\ 2 L2 (t)/dt = —2 k \\dp/dx\\ 2 L2 (t), characterizing Fickian diffusion. 

Next, consider a bounded domain V = [0, L] and let us suppose (as in 
the previous Sections) that at x = 0 zero-flux conditions are enforced. Eq. 
(1381) now becomes 

J t (\\p\\h(t) + j || J|| U*)) = -4o».tey - l IIJHi. W . (41) 

where the boundary term /boundary equals 

/boundary = 2 p(//, t) e/(//, t) 

= 2c [p + {L,t) +p-(L,t)\ [p + {L,t)-p~(L,t)\ . (42) 

Irreversibility, or equivalently the dissipative nature of the process, imply 
that the r.h.s. of eq. (1411) should be strictly non-positive. Consider at x = L 
homogeneous boundary conditions for the partial waves, such as in eq. (USD. 
A sufficient (not necessary) condition is that the constant 7 entering eq. (1301) 
satisfies the inequality 7 < 1, which is the condition obtained in Section [5] 
from positivity arguments. 

6 Biased hyperbolic transport and boundary-layer 
polarization 

In this Section we extend the analysis of boundary conditions in hyperbolic 
models to biased transport. Specifically, let v(x) be a deterministic velocity 
field, and consider the stochastic dynamics: 

dx(t) = v(x(t))dt + c (—l) x ^ dt , (43) 
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defined in the interval x € [0, L]. The PPM that specifies the statistical 
properties of eq. (|43)) is expressed by the system of two balance equations 


dp + (x, t ) 
dt 

dp~{x,t ) 
dt 


d[v(x)p + (x,t)j dp + (x,t ) 

dx dx 

d[v(x)p~(x,t)] dp~(x,t) 

dx +C dx 


— ap + (x, t) + ap (x,t) 
+ ap + (x, t ) — ap~(x, t ) 


(44) 


where p ± {x,t) are defined as in eq. (1491) . 

Let us study the evolution of the overall concentration field p(x, t ) and of 
its partial waves in a closed system, supposing that u(0) = 0, while v(L) > 0. 
For instance, consider the simple model for the biasing field 

V ( X ) = - y L sin > ( 45 ) 

where Vl > 0. Since v(L) > 0, there is a nonvanishing, outwardly directed, 
convective flux at x = L. The closedness of the system implies that the net 
flux vanishes at x = 0, L. Since u(0) = 0, reflecting boundary conditions for 
the partial waves apply at x = 0, i.e. eq. as in the previous examples. 

More interesting is the analysis at x = L, where the vanishing of the net 
flux dictates 

v(L)p(L,t) + J(L, t) = 0 , (46) 

where J(L , t) is the “diffusive” (non-Fickian) flux associated with the stochas¬ 
tic Poisson perturbation in eq. (1431) . In terms of partial waves this implies 

v{L) [p + (L,t) +p~(L,y)\ T c [p + {L,t) - p~{L,y)\ = 0 , (47) 

meaning that the “diffusive” flux, proportional to p + — p~, should counter¬ 
balance the convective one, proportional to p + +p~, creating a concentration 
boundary layer ( polarization ) near x = L. Since v(L) = Vl, the latter con¬ 
dition becomes 

(c - V L ) p~ (L,t) = (c + V L )p + {L, t) . (48) 

Eq. (1451) applies, without returning negative values of p~(L,t), if the con¬ 
dition 

c>V L , (49) 

is fulfilled. The case c < Vl provides negative values of p~(L,t), while if 
c^-Vl the backward wave at the boundary diverges, p~(L,t) oo. 

Therefore, boundary layer polarization can be properly defined in hy¬ 
perbolic transport models if eq. (1491) is fulfilled, i.e., the intensity of the 
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stochastic velocity perturbation, namely c, is strictly greater than the ve¬ 
locity V L at the boundary. In this case, the proper boundary condition for 
the partial waves at x = L is homogeneous and equal to 

p-(L,t) = Tp+(L,t) , T = ^±, 1 < T < oo . (50) 

c-V L 

Let us consider a numerical example which highlights better the meaning 
and the nature of condition (H9l) . We assume for v{x) the expression (H5ll 
and let L = 1. Figure [7] (panel a) depicts the behavior of the steady-state 
distribution p s (x) = pf(x) + pj(x) associated with eq. (HU), obtained from 
the asymptotic solution of PPM eq. (HU for t —>• oo. The values of the 
parameters are r c = 0.5, and k = 1, and we set Vl = vc with v < 1. The 
steady-state distributions refer to several values of v = Vl/c. Given r c and 
k, the parameters c, a specifying the stochastic perturbation are determined 
by the relations a = 1/2r c , c = y/2ka = \fkjr c . The numerical results of 
the PPM eq. (1441) . represented by solid lines in figure [71 are compared with 
stochastic simulations by considering an ensemble of N = 2 x 10 5 particles 
following the stochastic micro-dynamics (1431) . The initial positions of the 
particles are uniformly distributed in the interval [0,1]. Reflective boundary 
conditions are assumed at x = 0 and x = 1. Whenever a collision at the 
boundaries x = 0 or x = 1 occurs, say at t = t*, the stochastic process is 
updated. This means that at t = t*, the factor (—l)*w changes sign, namely, 
(_l)xh+) = — (—l) x (*-), where t± = lim^oi* ± e. As can be observed 
the agreement between PPM and stochastic simulations is excellent. The 
hyperbolic transport model provides the occurrence of a polarization layer 
near x = 1, that becomes more pronounced as v tends to 1. 

The maximum of the steady-state distribution occurs at some interme¬ 
diate abscissa x*, corresponding to the stable node of the forward wave, i.e., 
v(x*) + c = 0, and not at x = 0 which represents the stable equilibrium 
point of the deterministic dynamics. 

Keeping track of the values of (—1)*^), it is possible to extract from 
stochastic simulation data the distribution of the two-partial waves pf(x) 
and p~(x), using their definitions (fl8l) as joint probabilities with respect to 
the dichotomous outcome of (— 1) X W The comparison between PPM and 
stochastic simulations for the partial waves is depicted in figure [7| (panel 
b), revealing also in this case an excellent agreement, limited solely by the 
relatively small size, N = 2 x 10 5 , of the particle ensemble considered in the 
statistics. 

The comparison of PPM simulation and stochastic data reveals an in¬ 
teresting property of the condition (1491) . In point of fact, there are neither 
physical limitations nor fundamental consistency constraints in simulating 
an ensemble of particles moving according to eq. (|43l) beyond the limit im¬ 
posed by eq. m, i-e., for Vl > c, imposing reflective conditions at the 
boundaries. 
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Figure 7: Panel (a): Unique-steady state probability density function p s (x) 
vs x for t c = 0.5, k = 1, at different values of u. Solid lines are the results 
of the partial wave model (1441) . symbols represent stochastic simulations for 
an ensemble of particles evolving according to eq. (1431) . Line (a) and (□): 
v = 0.8, line (b) and (■): v = 0.9, line (c) and (•): v = 0.95. Panel (b): 
Steady-state probability density function p s (x) vs x (line (a) and symbols 
(□)), and associated steady-state partial waves pf{x) (line (b) and symbols 
(■)) and Px(x) (line (c) and symbols (•)) for r c = 0.5, k = 1, v = 0.95. 


While eq. (1281) stems from physical principles (non-negativity of proba¬ 
bilities, or equivalently dissipativity as developed in Section [5]), this is not 
the case of eq. (1491) . The origin of the critical condition c = Vl is not 
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of basic physical nature, but possesses a dynamic explanation. Namely, it 
is essentially a bifurcation point of the partial-wave model, associated with 
the ergodicity-breaking [28], corresponding to the birth of multiple invariant 
measures. For c > Vl there is a unique invariant stationary probability den¬ 
sity function p s (x). Conversely, for c < Vl, the stochastic dynamics (1431) is 
characterized by the occurrence of two disjoint invariant intervals I\ = [0, x*] 
and I 2 = [x**, 1], where x* is the first root of the equation v(x*) + c = 0 and 
x** the root of v(x**) — c = 0 (see figure [8]). Consequently there exist two 
distinct steady-state distributions p s ,h( x ), h = 1,2 localized within the two 
intervals //, . 



Figure 8: v(x) + c (line a) and v(x) — c (line b) for v = Vl/c = 1.1. Dotted 
vertical lines mark the two characteristic points x*, x** representing the 
endpoints of the two invariant intervals I± = [0, x*] and I 2 = [x**, 1]. 

This phenomenon can be easily verified via stochastic simulations start¬ 
ing from different initial distributions, as depicted in figure [9] panels (a)-(c) 
for v = 1.1. In this case, x* ~ 0.2412 and x** — 0.9087. 

Panel (b) depicts the steady-state distributions obtained by considering 
an initial particle ensemble uniformly distributed in [0,0.3] C I±, returning 
an absolutely continuous invariant distribution localized near x*. Conversely 
for an initial ensemble contained in I 2 (depicted in panel (c)), an impulsive 
invariant distribution p s , 2 (%) localized at x = 1 is obtained. The simulation 
starting from an uniform distribution throughout [0,1], depicted in panel (a), 
provides a mixed steady-state distribution p s {x) = ap s ,i(x) + (1 — a)p s ^{x), 
0 < a < 1, that is a convex combination of the two p s ,h( x )- 

Although a thorough analysis of ergodicity-breaking bifurcations is de¬ 
veloped elsewhere [28], since it goes beyond the analysis of boundary condi- 
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tions representing the focus of the present work, it is interesting to analyze 
the boundary conditions that apply to the partial-wave model for Vl > c in 
order to reproduce the results obtained from stochastic simulations. Also in 
the analysis of this issue, the answer stems from a wave-oriented interpreta¬ 
tion of hyperbolic transport. 

With the aid of figure[8l a wave-dynamic interpretation of the bifurcation 
occurring at Vl = c can be achieved. From their definition and dynamics 
dMD, the partial probability p + (x,t ) and p (x,t) corresponds to waves that 
propagate with velocities u + (x) = v{x)+c and v~(x) = v(x)—c, respectively. 
For Vl < c, in the neighbourhood of x = 1, v + {x) > 0 and v~(x) < 0, thus 
representing properly a progressive and a regressive wave, respectively. 

In this case, the balance at x = 1 of the overall flux permits to derive 
the condition (1501) for the “reflection coefficient” T, providing a regressive 
wave that balances exactly the outward probability flux carried out by the 
progressive wave. 

The bifurcation point Vl = c marks the transition for p~(x,t ) from a 
regressive wave ( v~(x ) < 0) to a progressive wave (u _ (x) > 0), near x = 1. 
In the latter case, impermeability at x = 1 cannot be longer enforced as 
relation connecting p~(l,t) to p + (l,t), since p~(l,t) would be negative. 
It is rather a consequence of the implicit assumption assumed in stochastic 
simulations of an arbitrarily large (and in principle infinite) potential barrier 
at x = 1 originating the reflective condition for eq. (1431) at x = 1. 

In order to derive the proper boundary conditions it is therefore con¬ 
venient to consider the restriction of eq. (1441) to each of the two invariant 
intervals Ih, h = 1,2. Consider eq. (jM]) restricted to x € I\ = [0, x*]. At 
x = 0, the usual reflective condition applies. At the other endpoint x = x* 
of the invariant interval, zero net-flux condition should be enforced, namely, 

[u(x*) + c] p + (x*,t) + (u(x*) — c] p~(x*,t) = 0 . (51) 

Since v(x*) + c = 0, eq. (1511) implies 

P~(x*,t) = 0, (52) 

which indicates that the regressive wave should vanish at x*. Figure [TO] de¬ 
picts the comparison of the PPM (l43l) restricted to Ii (solid lines) equipped 
with the boundary condition (1521) and the results of stochastic simulation 
for p S) i(x) and its partial components p+^x), p~ 1 (x) at v = 1.1. An en¬ 
semble of = 10 6 particles has been used in the stochastic simulations. As 
expected from eq. (ESI), stochastic simulation data confirm the vanishing 
value of the regressive wave at x*. 

Consider xE/ s = [x**, 1]. Similarly to eq. (1521) one derives 

p + (x**,t) = 0, (53) 

since v~(x**) = v(x**) — c = 0. 
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Figure 9: Steady-state distributions for r c = 0.5, k = 1, v = 1.1 (i.e., vi¬ 
olating condition (I49p f obtained from stochastic simulation of eq. (1431) for 
different distributions of initial conditions. Panel (a): uniform distribution 
in [0,1], panel (b): uniform distribution in [0,0.3], panel (c): uniform dis¬ 
tribution in [0.92,1], 


As stated above, the boundary condition at x = 1 does not derive from 
a balance between partial waves (as both are progressive waves in the neigh¬ 
bourhood of x = 1), but is a consequence of the assumption of an infi¬ 
nite potential barrier preventing particles to escape from the boundary in 
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Figure 10: Steady-state distributions p Si i(x), pf^x) restricted to I\ for r c = 
0.5, k = 1, v = 1,1. Solid lines represent the result of the partial wave model 
enforcing condition (15211 , symbols correspond the stochastic simulation data. 
Line (a) and (□): p St \(x), line (b) and (■): p+^x), line (c) and (•): p“ 1 (x). 


stochastic simulations. 

This observation implies that the partial waves possess an absolutely 
continuous component p±(x,t) and an impulsive part localized at x = 1, 

p ± (x, t ) = p^ix, t ) + ^(t) 5(x - 1) , (54) 

where the values of ^(t) stem from probability conservation. Asymptot¬ 
ically, a steady impulsive distribution is obtained, p S) 2 ( 2 ) = H x ~ 1)> as 
physically expected since v ± (x) > 0 uniformly in I 2 , as can be observed in 
figure [9] in panels (a) and (c). 

7 Concluding remarks 

In this article we have analyzed in detail the structure and the admissibility 
of different forms of boundary conditions for one-dimensional hyperbolic 
transport problems on the interval, both in the absence and in the presence 
of a deterministic biasing field. 

We stress again that the concept of a hyperbolic transport model should 
be regarded as a system of balance equations, obtained from the underlying 
micro-dynamics expressed in the form of a stochastic differential equation 
driven by dichotomous noise and its generalizations. In this framework, hy¬ 
perbolic transport theory is the wave-like counterpart of the classical trans- 
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port approach in which the stochastic micro-dynamics is controlled by a 
Langevin equation driven by Wiener processes (Brownian motion). 

In this concluding Section we would like to pinpoint some general ob¬ 
servations emerging from the present work that can be used to develop a 
self-consistent theory of hyperbolic transport, not limited to one-dimensional 
problems, and that be will thoroughly elaborated elsewhere: 

1. The primitive observables in hyperbolic transport theory are the par¬ 
tial probabilities (in an one-dimensional model with a single velocity c, 
these are just the partial probability waves p + (x , t) and p~(x, t)). The 
existence of these observables derives intrinsically from the dichoto¬ 
mous nature of the underlying stochastic fluctuations at microscale. 

2. The overall concentration p(x , t ) as well as its flux J(x, t ) should be re¬ 
garded in the theory as derived quantities, although they represent the 
physical observables commonly measured in transport experiments. 
This result is somehow similar to the setting of non-relativistic quan¬ 
tum mechanics where the wave function is the primitive quantity of 
the theory, while the probabilistic description of a quantum system is 
associated to the square of its modulus (Born condition). 

3. Boundary conditions in the hyperbolic theory are written as relations 
connecting forward and backward probability waves, corresponding 
to transmission/reflection conditions for the partial probability waves. 
The set of admissible boundary conditions is dictated by the positivity 
of partial probability functions, and can be derived alternatively by 
invoking the dissipative nature of the dynamics (Section [5]). 

4. There is a close connection between the functional relation describ¬ 
ing dissipation in hyperbolic “diffusion” eq. (I4U1) and the approach 
followed in Extended Irreversible Thermodynamics dm nn H21, where 
the entropy function depends also on the flux and not only on the 
concentration. The observation that the basic primitive quantities in 
hyperbolic transport theory are the partial waves suggests that the 
thermodynamic formalism developed in the theory of Extended Irre¬ 
versible Thermodynamics could be reformulated and/or generalized in 
order to include the primitive variables of hyperbolic transport theory 
(namely the partial probabilities) as the system variables with respect 
to which the state functions (such as entropy) are expressed. This 
effort would permit to resolve the problem of the absence of any prob¬ 
abilistic meaning of the Cattaneo equation (which is the prototype of 
generalized transport equation with memory in extended thermody¬ 
namics) in space dimensions greater than one. 
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